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ABSTRACT 

The  motion,  velocity,  and  acceleration  responses  at  ship 
center  of  gravity  (CG)  in  long-crested  seas  were  predicted 
for  a preliminary  conventional  hull  design  of  a new  class  of 
nuclear-powered,  guided-miss i le  strike  cruiser.  Time  histories 
from  the  computer  predictions  are  stored  on  computer  permanent 
files  for  easy  access  and  future  use.  Recently  developed 
techniques  for  determining  ship  response  time  histories  at 
various  locations  other  than  the  ship  CG  and  In  short-crested 
seas  were  used  to  calculate  sample  time  histories  of  important 
ship  responses. 

A computer  program  was  developed  to  read  the  predicted  time 
histories  and  compute  probabi 1 i t ies  that  specified  motion  response 
limits  would  be  exceeded  for  various  time  periods. 

ADMINISTRATIVE  INFORMATION 

This  work  was  funded  by  the  Naval  Ship  Engineering  Center  under  Work 
Requests  WR  66038  and  WR  66039  dated  2k  May  1976  and  WR  75239  dated  10  November 
1976.  The  work  is  identified  at  the  David  W.  Taylor  Naval  Ship  Research  and 
Development  Center  (DTNSRDC)  as  Work  Units  1-1568-863  and  1 - 1 568- 866 . 

INTRODUCTION 

As  part  of  the  process  of  designing  a strike  cruiser  (CSGN)  class,  naval 
ordnance  designers  require  estimates  of  the  motion,  velocity,  and  acceleration 
of  the  ship  at  positions  where  guns  and  radar  control  antennas  will  be  located. 

A rather  complete  prediction  of  significant  ship  motions,  velocities, 
and  acceleration  was  made  for  two  preliminary  designs  of  theCSGN  in  various 
sea  conditions.  However,  the  predictions  necessary  for  ordnance  performance 
evaluation  require  actual  time  histories  of  the  ship  responses  at  antenna 
and  weapon  locations  from  which  statistics  about  the  duration  of  extreme 
motion  responses  can  be  determined.  Prediction  techniques  and  computer 
programs  do  exist  to  predict  ship  response  time  histories  at  the  ship  center 


of  gravity  (CG).  Programs  are  also  available  to  predict  response  time  histo- 
ries for  short-crested  seas  and  to  predict  response  time  histories  for  any 
location  on  the  ship  from  the  time  histories  at  the  ship  CG. 

In  the  present  investigation,  a short  supplemental  program  was  written 
to  analyze  any  given  time  history  of  response  and  to  tabulate  statistics  about 
the  amount  of  time  the  response  exceeds  a specified  amplitude  limit.  These 
programs  have  been  used  to  predict  ship  response  time  histories  for  the  con- 
ventional CSGN  hull,  and  samples  of  the  results  are  shown.  Copies  of  the 
programs  used  are  included  with  the  time  histories  on  the  computer  disk  packs 
which  are  on  file. 


DESCRIPTION  OF  HULL  FORM 

This  hull  configuration,  C 17.2,  is  exactly  the  same  as  the  conventional 
hull  used  In  the  previous  CSGN  seaworthiness  prediction.  Table  1 shows  the 
principal  dimensions  of  the  hull  as  calculated  from  the  ship  offsets  and  used 
by  the  computer.  Figure  1 is  a body  plan  of  the  ship  as  used  by  the  computer. 

A preliminary  set  of  bilge  keels  (not  shown  in  the  computer  body  plan)  was 
Included  in  the  mathematical  model  of  the  ship. 

The  computer-determined  dimensions  and  form  coefficients  agreed  well  with 
the  design  values.  For  example,  the  calculated  displacement  of  16,866  long 
tons  (17.136  tonnes),  not  including  bilge  keels  and  skeg,  was  only  about  1.8 
percent  less  than  the  design  displacement  of  17,172  long  tons  (17, 448  tonnes). 
Addition  of  the  displacement  of  bilge  keels  and  skegs  would  make  the  difference 
even  less. 


PREDICTION  PROCEDURE 


Ship  responses  were  computed  by  using  the  CDC  6700  and  the  CDC  6600 
computer  systems  available  at  the  Center.  All  of  the  ship  motion  predictions 
are  stored  on  two  844  removable  disk  packs  titled  CSGN01  and  CSGN02. 


*Motter,  L.E.  and  T.R.  Applebee,  "Seaworthiness  Predictions  for  Two  Preliminary 
CSGN  Designs,"  DTNSRDC  Ship  Performance  Department  Report  SPD  724-01  (Sep  1976). 
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Many  computer  programs  were  used  to  predict  the  ship  rigid  body  responses 
and  statistical  predictions  presented  in  this  report.  Most  of  the  major 
programs,  as  well  as  some  of  the  smaller  access  programs,  are  documented  and 
a reference  for  each  is  given  herein  as  it  is  discussed.  A program  necessary 
to  analyze  the  time  history  data  is  documented  in  the  appendix.  A copy  of 
this  analysis  program  is  included  on  the  disk  pack,  and  instructions  are 
given  in  the  appendix  for  accessing  the  disk  pack  stored  data  and  program. 

PRESENT  OPTIONS 

Figure  2 is  a flow  chart  of  the  present  options  for  predicting  time  histo- 
ries and  the  procedure  for  generating  each  type.  The  regular  wave  ship  motion 
response  amplitude  operators  (RAO's)  referred  to  in  the  top  of  the  chart  were 

generated  and  stored  on  a disk  pack  by  using  the  DTNSRDC  Ship-Motion  and  Sea- 

2 

Load  Computer  Program.  As  shown  in  Figure  2,  ship  motion  time  histories  in 
any  specified  long-  or  short-crested  sea  can  be  predicted  at  the  ship  CG  or 
at  any  other  location  on  the  ship  from  the  RAO's.  The  cosine  squared  spreading 
function,  as  recommended  in  Reference  3,  is  used  to  predict  ship  motion  in 
short-crested  seas.  Any  of  these  time  histories  can  be  generated  by  using  a 
program  called  Time  History  Access  Program  which  is  documented  in  Reference  4. 

Once  any  time  history  has  been  generated  and  stored  on  the  disk  pack,  a 
time  history  analysis  program  (limit)  is  available  to  access  it  and  tabulate 
the  number  of  times  a motion  time  history  exceeds  a preassigned  motion  ampli- 
tude limit  (referred  to  as  an  event),  to  determine  for  how  long  the  limit  is 
exceeded  for  each  event,  and  to  estimate  probabilities  that  the  limit  will  be 
exceeded.  A more  detailed  explanation  of  the  program  operation  is  given  in 
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Meyers,  W.G.  et  al.,  "Manual  - NSRDC  Ship-Motion  and  Sea-Load  Computer  Program," 
NSRDC  Report  3376  (Feb  1975). 

^Baitis,  A.E.  et  al.,  "Design  Acceleration  and  Ship  Motions  for  LNG  Cargo 
Tanks,"  presented  at  10th  Naval  Hydrodynamics  Symposium,  Mass.  Inst.  Tech. 

(Jun  1974). 
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Baitls,  A.E.  et  al.,  "A  Non-Aviation  Ship  Motion  Base  for  the  DD  963,  CG  26, 

FF  1052,  FFG  7,  and  the  FF  1040  Ship  Classes,"  DTNSRDC  Ship  Performance 
Department  Report  SPD-738-01  (Dec  1976). 
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the  appendix  Sample  results  of  the  time  history  analysis  are  presented  later 
In  this  report. 


TIME  HISTORY  PREDICTION 


i 


The  ship  motion,  velocity,  and  acceleration  RAO's  presented  in  Reference  1 
and  stored  on  disk  pack  CSGN01  are  used  as  input  to  the  time  history  prediction 
program.  Because  of  storage  requirements  involved,  only  time  histories  for 
those  conditions  listed  in  Table  2 were  generated  and  stored  on  disk  packs 
CSGN01  and  CSGN02 . 

For  each  condition  in  Table  2,  time  histories  for  pitch,  heave,  roll,  and 
vertical  acceleration  at  the  ship  CG  in  long-crested  seas  were  generated  as 
examples  of  the  normal  time  history  prediction. 

Bretschneider"’  Irregular  sea  spectra  o.  1.0-ft  significant  wave  height 
and  modal  wave  periods  of  9,  11,  13,  and  15  sec  were  used  in  the  prediction. 
Since  ship  responses  are  regarded  as  being  linear  with  wave  height  in  this 
work,  the  responses  in  seas  with  greater  wave  heights  can  be  obtained  by  taking 
the  results  for  unit  significant  wave  height  and  multiplying  by  the  desired 
significant  wave  height  In  feet. 

The  modal  wave  period  is  the  period  at  which  the  maximum  energy  or  peak 
of  the  wave  energy  spectrum  occurs.  The  time  history  predictions  cover  a 
<-ange  of  realistic  modal  periods,  including  some  most  likely  to  occur  at  sea 
as  well  as  some  that  are  less  likely  to  occur  but  stimulate  much  more  severe 
ship  motion  responses. 

Considerably  more  information  was  generated  by  the  computer  than  is 
reasonable  or  necessary  to  be  shown  in  this  report.  Accordingly,  only  those 
conditions  and  ship  responses  of  particular  significance  will  be  shown.  These 
results  will  serve  as  examples  of  the  available  methods  and  possible  applica- 
tions of  the  prediction  techniques. 

Figure  3 shows  an  example  of  a roll  and  a heave  acceleration  time  history 
as  generated  by  the  computer.  These  time  histories  are  for  a 30-knot  ship 
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Bretschne ider , C.L.,  "Wave  Generation  by  Wind,  Deep  and  Shallow  Water,"  in 
"Estuary  and  Coastline  Hydrodynamics,"  edited  by  Arthur  T.  Ippen,  McGraw-Hill 
Inc  , New  York  (1966)  pp.  I 33~ 1 96 . 
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speed  at  75“deg  heading  in  a long-crested  sea  with  1 6 . 9~ f t (5.15  m)  significant 
wave  height  and  a spectra  modal  period  of  9 sec.  This  condition  causes  one 
of  the  most  severe  responses  for  roll. 

To  illustrate  use  of  the  origin  shift  option  in  the  Access  Program,  time 
histories  were  generated  in  long-crested  seas  for  vertical  displacement, 
velocity,  and  acceleration  at  the  center  of  the  roller  path  plane  of  the  bow 
gun.  This  gun  will  be  located  on  the  ship  longitudinal  centerline  95.0  ft 
(28.96  m)  aft  of  the  forward  perpendicular  and  53.0  ft  (16.15  m)  above  the 
ship  baseline.  Note  that  since  roll  and  pitch  are  the  same  everywhere  on  the 
ship,  roll  and  pitch  time  histories  are  the  same  at  the  gun  as  at  the  CG. 
Finally,  roll  and  vertical  accelerations  at  the  CG  in  short-crested  seas  were 
generated  as  examples  of  the  short-crested  time  history  program.  The  headings, 
ship  speeds,  and  modal  periods  of  the  sea  spectrum  used  for  each  response  are 
the  conditions  which  give  the  most  or  nearly  the  most  severe  response  in  a 
sea  with  a l6.9"ft  (5.15  m)  significant  height. 

TIME  HISTORY  ANALYSIS 

Time  histories  of  ship  responses  have  been  predicted  and  analyzed  for 
various  ship  location  and  sea  condition  to  illustrate  how  the  time  history 
prediction  and  analysis  programs  can  be  used. 

All  sample  calculations  were  made  for  a ship  speed  of  30  knots  since 
motions  were  more  severe  than  at  the  two  lower  speeds  considered.  All  sample 
calculations  considered  the  ship  as  operating  in  a Sea  State  5 with  a signifi- 
cant wave  height  of  16.9  ft  (5.15  m) . The  modal  period  of  the  wave  spectra  is 
given  with  each  example. 

Table  3 enumerates  the  ship  responses,  head!ngs,  and  spectral  modal 
periods  for  which  time  histories  were  developed.  In  addition  to  presenting 
rms  values  for  each  response,  the  table  presents  four  columns  of  numbers  which 
appear  in  the  form  0.7/393.  The  first  number,  or  numerator  of  the  fraction, 
is  the  amplitude  limit  used  in  the  analysis  of  the  ship  response  time  history 
to  determine  the  number  of  cycles  which  exceed  the  limit.  The  second  number, 
or  denominator  of  the  fraction  titled  the  number  of  events,  is  the  number  of 
cycles  which  exceed  the  specified  amplitude  limit  in  30  min  of  ship  operation. 


The  four  limits  were  chosen  in  an  effort  to  ensure  that  no  less  than  about 
10  percent  nor  more  than  about  80  percent  of  the  cycles  In  a given  time  history 
would  exceed  any  of  the  limits.  Generally  about  10  percent  of  the  total  number 
of  cycles  in  a 30-min  time  history  must  exceed  the  limit  in  order  to  have  a 
sufficient  number  of  events  to  generate  consistent  statistical  properties  about 
time  durations  when  limits  will  be  exceeded.  At  the  other  extreme,  it  is 
assumed  that  the  resulting  duration  time  statistics  will  be  out  of  the  range 
of  useful  information  when  80  to  90  percent  of  cycles  exceed  a given  amplitude 
1 imi t . 

The  probability  of  exceeding  a ship  response  limit  for  each  of  the  ship 
responses  in  Table  3 is  shown  in  Figures  4a  through  4i  in  the  same  order  as 
listed  in  the  table.  The  curves  in  the  figures  represent  the  probability  dis- 
tribution function  as  determined  from  analyzing  a 30-min  time  history  several 
times  with  a different  limit  for  each  analysis.  The  individual  points,  from 
each  analysis  of  the  time  history,  were  so  consistent  that  only  the  faired 
distribution  function  curve  is  shown  in  the  figures. 

It  is  known  that  the  amplitudes  of  ship  motion  time  histories  are  Rayleigh 
distributed  and  have  a distribution  function  of 


Prob  { r > R } = e 2 a2 
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where  R^  is  the  response  limit  amplitude  and  a is  the  root  mean  square  (rms). 
Using  this  formula  and  the  frequency  domain  results  in  Reference  1,  the 
probabilities  of  exceeding  each  of  the  four  limits  in  Table  3 were  computed 
and  shown  in  Figures  4a  through  4i  with  the  time  history  analysis  results.  As 
seen  in  the  figures,  there  was  good  agreement  between  the  probabilities  obtained 
from  the  time  domain  analysis  and  from  the  frequency  domain  results.  Thts  would 
be  expected  since  the  time  histories  were  derived  from  the  frequency  domain 
results  and  is  mentioned  as  evidence  that  the  time  history  prediction  from 
frequency  domain  results  does  work  well.  The  figures  indicate,  for  example, 
that  If  a gun  cannot  operate  effectively  when  the  roll  is  greater  than  +_  16  deg, 
then,  as  shown  in  Figure  4c,  its  operation  will  be  interrupted  by  30  percent 
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of  the  waves  encountered  while  operating  at  30  knots  in  a 75"deg  heading,  and 
16.9  ft  (5-15  m)  significant  wave  height. 

The  percentage  of  the  total  run  time  that  the  response  exceeds  any  limit 
is  shown  in  Figures  5a  and  5i  in  the  same  order  as  listed  in  Table  3.  The  time 
was  computed  on  the  conservative  side  from  the  time  history  analysis  by  assuming 
that  each  event  of  the  0.5-sec,  time-duration  interval  occurred  for  the  upper 
time  limit  of  the  interval.  These  figures  indicate,  for  example,  that  the 
roll  angle  will  be  greater  than  +_  16  deg  for  about  12.5  percent  of  the  total 
run  time  even  though  about  30  percent  of  the  roll  cycles  will  exceed  +_  16  deg 
while  the  ship  is  operating  at  30  knots  and  a 75"deg  heading  with  respect  to 
the  waves  in  a State  6 sea. 

Figures  6a  through  6i  show  the  conditional  probability  density  functions 
of  the  time  durations  given  that  the  specified  limit  has  been  exceeded.  These 
functions  have  been  normalized  so  that  the  area  in  each  bar  of  the  graph 
represents  the  probability  that,  the  limit  having  been  exceeded,  its  duration 
is  within  the  duration  period  of  the  time  interval  of  the  bar.  Also,  the 
total  area  of  all  the  bars  in  each  graph  is  1.0,  i.e.,  given  that  the  response 
limit  was  exceeded,  there  is  a 100-percent  chance  that  it  was  exceeded  for 
some  length  of  time.  As  seen  in  the  figures  for  all  ship  responses,  as  the 
response  limit  is  increased,  the  most  frequently  occurring  time  duration  is 
shorter . 

As  an  example  of  the  meaning  of  these  curves, in  Figure  6c  the  probability 
that  the  roll  cycles  will  exceed  +_  17  deg  for  a duration  time  between  1.5  and 
2.0  sec  is  48  percent  per  second  times  0.5  sec,  or  24  percent  of  the  cycles. 

Figures  7a  through  7i  show  joint  probability  density  functions  of  the 
time  durations  for  four  different  amplitude  limits.  Again  the  density  functions 
have  been  normalized  so  that  the  area  in  each  bar  of  a graph  represents  the 
joint  probability  that  (l)  the  limit  is  exceeded  and  (2)  the  period  of  the 
exceedance  is  within  the  duration  time  of  the  time  interval  of  the  bar.  The 
total  area  in  all  the  bars  for  any  one  response  limit  is  the  probability  that 
the  response  limit  is  exceeded,  as  shown  in  Figure  4.  For  example,  consid- 
ering Figure  7c,  the  probability  that  a roll  cycle  will  exceed  +_  17  deg  for 
a duration  time  of  1.5  to  2.0  sec  is  12.5  percent  per  second  times  0.5  sec  or 
6.3  percent. 
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To  smooth  out  the  bar  graphs  and  to  enable  a comparison  of  the  joint 
probability  densities  shown  in  Figures  7a  through  7 • » each  bar  graph  in  Figure 
7 was  hand  faired  to  a smooth  curve.  An  attempt  was  made  to  maintain  the  same 
area  under  the  curve  as  in  the  correspondong  bar  graph.  The  resulting  faired 
curves  are  presented  in  Figures  8a  through  8i  to  show  the  relationship  between 
the  density  functions  for  different  response  limits.  As  seen  in  these  figures, 
the  curves  for  different  response  limits  vary  in  a regular  fashion  making  it 
possible  to  estimate  the  time  duration  probabi 1 1 1 ies  for  response  limits  other 
than  those  used  in  the  calculations. 

RECOMMENDATIONS 


5» 


P 
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The  cost,  time  and  computer  storage  required  to  predict  and  analyze  time 

histories  is  considerable.  If  the  mathematical  form  of  the  probability  density 

function  for  the  duration  time  could  be  developed,  it  may  be  possible  to  make 

predictions  of  the  duration  time  interval  directly  from  the  rms  results.  Work 

6 7 8 

has  been  done  by  Rice,  Price  and  Bishop'  and  by  Tikhonov  with  some  success 
on  simplified  cases  using  numerical  approximations.  With  the  examples  pre- 
sented in  this  report  as  a guide,  a reasonable  method  of  predicting  the 
duration  time  interval  from  the  rms  results  may  be  possible,  and  should  be 
attempted. 


Rice,  S.O.,  "Mathematical  Analysis  of  Random  Noises,"  N.  Wax,  Ed.  of  Selected 
Papers  on  Noise  and  Stochastic  Processes,  Dover  Publications,  Inc.  (1951*). 

^Price,  W.G.  and  R.E.D.  Bishop,  "Probabilistic  Theory  of  Ship  Dynamics," 
Chapman  and  Hall  ( 1 97^ ) • 

8 

Tikhonov,  V. I . , "The  Distribution  of  the  Duration  of  Excursions  of  Normal 
Fluctuations.  Non-Linear  Transformations  of  Stochastic  Processes,"  Ed., 
Kusnetsov,  Stratonovick  and  Tikhonov. 
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APPENDIX 
PROGRAM  LIMIT 


The  LIMIT  program  is  used  to  analyze  time  histories  that  have  been 
generated  using  the  Time  History  Access  Program  (THACP) , and  to  obtain  various 
probabilities  with  relatton  to  an  imposed  limit.  The  limit,  which  is  an  input 
variable  is  a boundary  applied  to  the  positive  and  negative  sides  of  a time 
history.  If  a cycle  exceeds  the  limit,  the  duration  of  time  above  the  limit 
is  noted.  The  predicted  time  histories  are  in  a digital  form  with  amplitude 
samples  every  0.5  sec  for  the  30-min  run.  Therefore,  time  durations  can  be 
determined  only  to  within  +0.25  sec.  For  this  reason  the  analysis  of  the 
time  histories  is  set  up  in  tabular  form  as  shown  in  Table  A.l.  The  number 
of  occurrences  of  exceedance,  referred  to  as  "events,"  are  ordered  according 
to  time  duration  periods  (see  example  of  output,  as  shown  in  the  first  two 
columns  of  Table  A.l).  From  this  information,  the  following  statistical 
percentages  are  calculated  and  shown  in  the  last  three  columns  of  Table  A.l: 

THE  CONDITIONAL  PROBABILITY  OF  DURATION  TIME,  defined  as  the  probability 
that  a particular  duration  time  will  occur  if  the  limit  is  exceeded 
(the  number  of  events/total  number  of  events). 

THE  TOTAL  TIME  EXCEEDED,  defined  as  that  part  of  the  entire  run  time 
that  the  limit  is  exceeded  for  a particular  duration  time  (the 
number  of  events  x time  duration/total  run  time). 

THE  JOINT  PROBABILITY,  defined  as  the  probability  that  for  any 
cycle,  the  limit  will  be  exceeded  for  a specific  time  duration 
(the  number  of  events/total  number  of  cycles). 

INPUT 

Data  Card  Set  1 , one  card,  FORMAT  (15) 

This  card  set  contains  one  integer  variable. 

NRSP,  Columns  1-5,  Is  used  to  determine  the  number  of  responses  to  be 
analyzed  (i.e.,  the  number  of  cards  to  be  read  in  Data  Card  Set  2). 


Data  Card  Set  2,  NRSP  cards,  FORMAT  (I5.5X.F10. A) 

This  card  set  contains  one  Integer  and  one  floating  point  variable. 
IRESP,  Columns  1-5*  is  the  response  sequence  number  as  defined  In  the 


output  of  the  Time  History  Access  Program. 

LIM,  Columns  11-20,  is  the  boundary  or  limit,  In  appropriate  units  to  be 
Imposed  on  the  time  history. 

(Note  that  Data  Card  Set  2 Is  repeated  for  each  new  response  or  limit  to  be 
analyzed.) 

A sample  of  the  entire  computer  card  deck  necessary  to  access  the  disk 
pack  and  run  the  program  is  shown  in  Figure  A. 1. 
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Figure  6 - Conditional  Probability  Density  Functions  for  Tine  Duration 
Given  the  Response  Limit  Is  Exceeded  for  30~Kriot  Speed  in 
Sea  State  6 (16.9  ft  (5.15  m)  Sianificant  Wave  Height) 


: IMF  mHATies  is  sfconos 

Figure  6a  - Pitch  in  Long  Crested  Waves 

( 1 8 0 deg  Heading,  15  sec  Modal  Period) 


PFUCEST 


0.22  ft'f  MM  IT 


0. ?'  a * ! T w T T 


l 1 ill  l 

*01  7 T 

Tit*  miATinN  is  srcnsn* 


i 


Figure  6d  - Heave  Acceleration  In  Lonq  Crested  Waves 
(105  deg  Heading,  9 sec  Modal  Period) 
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Figure  7 " Joint  Probability  Density  Function  Bar  Graphs  for  Time 
Duration  and  Response  Limit  for  30~Knot  Speed  In  Sea 
State  6 (16.9  ft  (5.15  m)  Significant  Wave  Height) 
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OOOOOOOOOOOOOOOOOOOOOO  ENU  UF  RECORD 


C*  DEBUG 

C*  ARRAYS 

PROGRAM  LIMIT  (INPUT, OUTPUT, TAPtS* I NPUT , T APEfeaOuTPUT , T APE50 • 

2 OEBuGaOUTPUT ) 

COMMON  /ZC/  SAMPLE (3600  I .NCYCLS 

DIMENSION  THI  ST  (19,'.,  30)  .TIMES  (21)  . 1SIGNI3600)  .UNITS  (6) 

DIMENSION  FPROBI21) ,TPROB<20) ,NUCC(21) ,\TP(3600)  , PROBJ(2D 
DIMENSION  TITLE  (8)  .SIGBHIR)  .TMOOAL  (*. ) .CT I TL(  19) 

DIMENSION  RESPTh ( 50,30) . T I TL ( 8 ) .PNT I TL ( 3 .SO ) ,RST I TL ( 3 .SO ) , 

2 IPOINT(SO)  . IMOTN(SO) . I TYPE ( 50 ) . 1PNT ( 50 ) . X AP ( 50 ) .YCLI50) .ZBL(SO) . 

3 JBLBLI50) .XSTAR(SO) .YSTAR(SO) .ZSTAR(SO) ,XMU(S0) 

REAL  LIM.MEAN 

DATA  TIMES  /0.0.0.S,  1 .0. 1 .5 ,2 .0 . 2 .5 . 3 .0 , 3 .5.A .0 .9 .S.S.O ,5 .S .6 . 0 . 
2 6.5,7.0,7.5.8.0,8.5.9.0,9.5.10.0/ 

DATA  UNITS  / 1 OHOEGREES  . 

2 10HFEET  . 

2 10HOEG/SEC  . 

2 10MFT/SEC  , 

2 10HOEG/SEC»«2. 

2 I OHG • S / 

READ  (5. 1000)  NRSP 
•RITE  (6,2005) 

DO  500  I RE  = 1 .NRSP 

READ  (50)  UTLE. HEAD. SPEED. NSPEL.SlGBH.THODAL, PI, X2P1.RD.NB.DB, 

2 NCHN.NREC.KPTS.DT .SR.TRUN.NPTS.CTI TL , BVSLOP . TPST ,ELL 

READ  (50)  NR^SP.TITL.SIGBHT.PERMOO.NPUINT.IPNT.PNTITL.XAP.YCL. 

2 ZBL.DBLBL.DAPO.XSTAR.YSTAR.ZSTAR.RSTITL.IMUTN, I TYPE. I POINT, xmu 
JO  5 1*1,21 
5 NOCC(I)  a 0 

READ  (5.1050)  IRESP.LIM 
I UN  IT  = 2 

IF  ( ITYPE ( IRESP) .EO.l .AND. IMOTNl IRESP) .GT.3)  1UN1T=1 

IF  ( ITYPE ( IRESP) .EQ.’.ANO. IMOTNl IRESP) .GT. 9)  IUNITa3 
IF  ( I TYPE ( IRESP) .EO. 2. ANO. IMOTNl IRESP) .LE. 9)  I UN I T = 9 
IF  ( ITYPE ( IRESP) .EO.l. AND. IMOTNl IRESP) .GT. 15)  IUNIT=5 
IF  ( I TYPE ( IRESP) .EO. 3. AND. IMOTNl IRESP) .LE. 15)  IUNIT=6 
M ■ 0 
10CC  a 0 
NTOCC  a 0 
K ■ 0 

IMPROB  a 0.0 
TIME  S 0.0 
MEAN  a 0. 

DO  100  NRal.NREC 

IF  (KPT  S*NR.GT . NPTS ) KPTS=MOD (NPTS.KPTS) 

READ  (50)  ( (RESPTHI I , J) , 1*1 .NRESP) . J« 1 .KPTS) 

DO  50  N* 1 .KPTS 
K « K.  1 

sample  (K)  a RESPTHI IRESP. N) 

MEAN  a MEAN.SAMPLE <K) 

IF  (SAMPLE  IK) .GE. 0.0)  ISIGN(K)=1 
IF  (SAMPLE (K) ,LT. 0.0)  ISIGN(K)=2 
IF  (K.EQ.l)  GO  TO  31 


igure  A1  - Listing  of  Limit  Program  with  Control  Cards  and  Data 
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IF  (ABS(SAMPLC<R-1>  ) .LE.LlM.AND.ABStSAMRLE (K>  > .GT.L1M)  00  TO  15 
IF  (A8S(5AHPLE(A-ln  .GT.LIM.ANU.AbSISAMPLEIK)  ) .GT.LIMJ  GO  TO  20 
10  IF  (ABS1SAMPLE <*-l > ) .GT.LIM.AND.ABSISAMPLE IK) > .lE.LIM)  GO  TO  25 
GO  TO  50 
15  TIME  « 0.5 
lOCC  « 1 
00  TO  50 

20  IF  (ISIGN(K) .NE.ISIGN(K-l) ) GO  TO  25 
TIME  » TIME. 0.5 
IOCC  * IOCC.1, 

GO  TO  50 

25  IF  (TIME. GT. 10.0)  GO  TO  35 
00  30  1*1.20 

IF  (TIME. GT. TIMESlI ) .AND. TIME. LE.T 1ME5I 1.1 ))  NOCC ( I > »NOCC ( I ) » 1 

30  CONTINUE 

31  IOCC  ■ 0 
TIME  * 0.0 
00  TO  50 

35  NOCC (21)  = NOCC ( 21 ) ♦ 1 
M ■ M.  1 

ATP  (M)  = (TIME/TRUN1M00.0 
lOCC  * 0 
TIME  = 0.0 
50  CONTINUE 
100  CONTINUE 

MEAN  * MEAN/K 
CALL  CYCLES  IMEAN.K) 

00  110  J= 1 . 21 
NTOCC  = NTUCC.NOCC ( J) 

110  CONTINUE 

IF  (NTOCC. EQ.O)  GO  TO  145 
00  120  J=  I .21 

pmobj(J)  = (Float(nocc« j) ) /float (ncylls) >*ioo.o 

TMPROB  = TMPROB  ♦ PROrtJ(J) 

120  FPROH(J)  = (FLOAT (NOCC(J) I /FLOAT (NTOCC) ) •100.0 
00  130  L=2 • 2 1 

130  IPROH(L-l)  = ( T IMES  (L  ) •FLOaT  (NOCC  (L-I  ) I /TRIJN)  *1  00.0 
00  140  I * 1 .M 
140  2 TP  * 2TP.XTP ( I ) 

145  WRITE  (5.1500)  T 1 TLE .HE AD . SPEED . S I OWHT . PE RMOO 

■RITE  (5.1500)  (RSTITL ( J.IRESP) . J=1 .3) .LIM. UNITS! IUNIT) 

IF  (NPOINT.EQ.O)  GO  TO  148 
I PT  = IPOINT ( IRESP) 

IF  (NSPEC.LT. 0)  WRITE  (8.1540) 

IF  (NSPEC.GT.O)  WRITE  (6.1530) 

■RITE  (6.1650)  (PNTITL(I.IPT) ,1=1,3) 

GO  TO  14R 

148  IF  (NSPEC.LT. 0)  WRITE  (6,1520) 

IF  (NSPEC.GT.O)  WRITE  (6.1510) 

149  IF  (NTOCC. EQ.O)  GO  TO  160 
WRITE  (6.1700) 

00  150  1=1.20 

■ RITE  (6.1800)  TIMES! I ) .TIMES! I *1 ) .NUCCI I ) .FPROBI I ) .TPROBI I I 
•PROBJ(I) 

150  CONTINUE 

■RITE  (6.1900)  NOCC(2l) ,FPRUH<21| .ZIP,  PRuHj(21) 

160  IF  (NTOCC. EQ.O)  WRITE  (6.2010) 

■RITE  (6.2000)  TRUN.SR.NCYCLS.TMPR08 
REWIND  50 
500  CONTINUE 
looo  format  (I5> 

1050  format  i is,5x,fio.4) 

1500  FORMAT  (•1»,//.21*.8A10.//.35X.*MEAUING  = *.F5.1.*  OtGREES»»5X. 

2 ’SPEED  * *.F4,1,a  KNOTS* ,/. 3SX , #S 1G.  « M s *,F5.2.*  FEET*.8X. 

3 *MOD.  PER.  = • »F4 . 1 . * SECONDS*) 

1510  format  (*o*.//40x.*—  longcrestuj  time  history  statistics  — 

2 /60x.*AT*./56X.*TMF  ORIGIN*) 
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1520  FORMAT  (*0*.//40X.*~  SHORTCRESTED  TIME  HISTORY  STATISTICS  — 

2 /61x,*AT*»/57x.*THE  ORIGIN*) 

1530  FORMAT  <*0*.//40x.*~  LONGCRESTED  TIME  HlstoRY  STATISTICS  — 

2  /60X.*AT*> 

1540  FORMAT  (*0*.//40x.*--  SHORTCRES1EU  I I HE  HISTORY  STATISTICS  — •. 

2  /61X.*AT*) 

1600  Format  (*0*. 34X, ‘RESPONSE  = *3A 1 0 . /35X , *L IMI T «*F6. 2 • I X . Al 0 ) 

1650  Format  (•  *.42x.taio) 

1700  FORMAT (*0*.44X, ‘NUMBER*. 5X.*COND.  PRUB.  OF  * . 4 X . *ToT AL  TlME*.7x, 

. *joint*/23x,*time  Dorati  on*  * 1 1 x **of  *• 7x,*ourat I ON  TImE*.6X. 
.•EXCEEDE0*.5X. ‘PROBABILITY*/  25X ,*( SECONDS )*. 1 1 X . ‘EVENT S* . 7X . 

. •< PERCENT ) *.BX.* (PERCENT) *.SX,* (PERCENT)*/) 

1800  Format  (•  *.23X,F4.1,*  - *«F4.) ,UX,I3,10X,F6.2.10X,F6.2,8X.F6.2) 
1900  FORMAT  (•  * . 25X . • > 1 0 . 0* . 1 3x , I 3 . I OX ,F 6 .2 , 1 OX ,F 6. 2 . 8X .F 6 .2 ) 

2000  format  (*o*,19x,*note:  total  ron  time  is  *,F6.I. 

2 • SECONDS  at  A sample  RATE  of  *.F3.1.*  POINTS  PER  SECOND. *./ .27x . 

3 ‘TOTAL  NUMBER  OF  CYCLES  l S* • 15. * . * . / ,27 X ,«THE  LIMIT  IS  EXCEEDED  • 

4 *BY  *.F6.2.*  PERCENT  OF  TmE  CYCLES.*) 

200S  format  (*i*.sx.*def initions:  nomber  of  events  * the  number  of  • 

2 ‘TIMES  THE  LIMIT  AMPLITUDE  IS  EXCEEDED  FUR  THE  SPECIFIED  TIME* 

3 • duration.*  //20x,*conoitiunal  probability  of  duration  time  ■ 

4 I HE  PROBABILITY  THAT  ThE  SPECIFIED  DURATION  TIME  WILL  OCCUR.  • 

5 /63X . *G  l VEN  VAT  THE  LIMIT  AMPLITUDE  HAS  BEEN  EXCEE0ED.*//20X . 

6 ‘TOTAL  TIME  EXCEEDED  = THE  PERCENT  UF  THE  TOTAL  RUN  TIME  THAT  THE 

7 LIMIT  AMPLITUDE  IS  EXCEEDED*/41x*  FOR  THE  SPECIFIED  DURATION* 

R • T IMF ,*//20x .‘JOINT  PROBABILITY  * THE  PROBABILITY  THAT,  FOR  ANY* 

b • cycle,  the  limit  amplitude  will  be  exceeded  and*/40x, *the* 

9 • Duration  time  will  be  as  specified. *i 
2010  format  i *o*////.ssx . i ih**  sorry  •*,//. 4 ox,*no  cycles  were  found  • 

2 • TO  EXCEED  THE  LIMIT*//////////) 

STOP 

END 

SUBROUTINE  CYCLES  (MEAN.NP) 

REAL  MEAN 

Common  /zc/  sample ( 3600) .ncycls 
DATA  EPS  /l.E-10/ 

NCYCLS  = 0 
JO  SO  IP= 1 .NP 
IF  ( IP.GT.l ) GO  TO  10 
PEA*.  = SAMPLE  ( 1 ) 
trough  = pea* 

SGNPRE  = SIGN ( I . .SAMPLE ( I ) -MEAN) 

10  CORR  = SAMPLE ( IP) -MEAN 

IF  ( ABS (CURR) .LT .EPS)  GO  TO  50 
SGNCUR  = SIGN( 1 ..CURR) 

IF  (SGNCUR. NE. SGNPRE ) GO  TO  30 
IF  (CURB. LT. 0.1  GO  TO  20 
IF  (SAMPLE ( IP) .LT.PE»*)  GO  TO  50 
PEA*  * SAMPLE (IP) 

GO  TO  50 

20  IF  (SAMPLE ( IP) .GT. TROUGH)  GO  TO  SO 
TROUGH  r SAMPLE (IP) 

GO  TO  50 

30  SGNPRE  = SGNCUR 
NCYCLS  3 NCYClS-1 
50  CONTINUE 
RETURN 
END 

0000000000000000000000  END  OF  RECORD 
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OOOOOOOOt  .000000000000  END  OF  FILE 
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TABLE  1 - COMPUTER-PREDICTED  PARTICULARS  FOR  CONVENTIONAL  HULL  C 17-2 


Displacement,  S.W. , long  tons  (Tonnes) 

16,866 

(17,136) 

Length  between  Perpendiculars,  ft  (m) 

666.0 

(203.0) 

Beam,  ft  (m) 

76. A6 

(23-31) 

LCG,  aft  of  Midship,  ft  (m) 

8.A7 

(2.58) 

Transom  Width,  ft  (m) 

A!  .AO 

(12.62) 

Waterplane  Area,  ft  (m  ) 

37,279 

(3A63) 

Vertical  Center  of  Gravity,  KG,  ft  (m) 

28.00 

(8.53) 

Transverse  Metacenter,  GM,  ft  (m) 

7.38 

(2.25) 

Vertical  Center  of  Buoyancy,  KB,  ft  (m) 

13-01 

(3-97) 

Vertical  Metacentric  Height,  KM,  ft  (m) 

35.38 

(10.78) 

Block  Coefficient,  Cb 

0.52 

Midship  Section  Coefficient, 

0.89 

Prismatic  Coefficient,  C 

P 

0.59 

TABLE  2 - CALCULATION  CONDITIONS  (SPEEDS,  HEADINGS  AND  SEA  CONDITIONS) 


Ship  Speeds,  knots 


5 


20 


30 


Headings,  deg 


0 15  30  45 


60  75  90  105  120  135  150  165  180 


Wave  Spectra  Modal 

Perod,  sec  9 II  13  15 


NOTES:  (A)  All  wave  spectra  have  a 1.0-ft  (0.3048  m)  significant  wave 

height  (double  amplitude). 

(B)  Motion,  velocity  and  acceleration  rms  responses  can  be  obtained 
for  any  sea  state  by  multiplying  the  predicted  rms  values  by 
the  significant  wave  height  of  the  sea. 

(C)  180-deg  heading  is  head  seas. 
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TABLE  3 - RESPONSE  AMPLITUDE  LIMITS  AND  NUMBER  OF  CYCLES  WHICH  EXCEED  THE  LIMIT 


180-deq  neadlnq  is  head  seas 


TABLE  A. 1 - SAMPLE  OUTPUT  FOR  PROGRAM  LIMIT 
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total  DUMBER  of  CYCLES  IS  305 


DTNSRDC  ISSUES  THREE  TYPES  OF  REPORTS 


(1)  DTNSRDC  REPORTS,  A FORMAL  SERIES  PUBLISHING  INFORMATION  OF 
PERMANENT  TECHNICAL  VALUE,  DESIGNATED  BY  A SERIAL  REPORT  NUMBER 

(2)  DEPARTMENTAL  REPORTS.  A SEMIFORMAL  SERIES  RECORDING  INFORMA 
TION  OF  A PRELIMINARY  OR  TEMPORARY  NATURE,  OR  OF  LIMITED  INTEREST  OR 
SIGNIFICANCE,  CARRYING  A DEPARTMENTAL  ALPHANUMERIC  IDENTIFICATION 

(3)  TECHNICAL  MEMORANDA,  AN  INFORMAL  SERIES.  USUALLY  INTERNAL 
WORKING  PAPERS  OR  DIRECT  REPORTS  TO  SPONSORS.  NUMBERED  AS  TM  SERIES 
REPORTS,  NOT  FOR  GENERAL  DISTRIBUTION 


